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1 Introduction 

Over the last years, highly sensitive magnetic field 
sensors, SQUIDs, have been employed to measure 
the magnetocardiogram (MCG), which are 
generated by ionic currents flowing in myocardium 
when it is active. 

The MCG method provides additional information 
in comparison to surface electric potential 
distribution. The obtained information is 
complementary in nature because the magnetic field 
is related to the curl of the impressed current sources 
into the human heart while the electric field is 
related to its divergence [1], It is known that the 
myocardium has many kinds of anisotropy. One of 
them is determined by an anisotropy of conductivity 
(uniform anisotropy), which plays a key role in 
realistic propagation of the electrical excitation in 
the myocardium. This was indicated by computer 
simulation of normal MCG based on a propagation 
model of an excitation wave front in three- 
dimensional arrays of cardiac cells within the 
cellular automata theory [2], Animal studies and 
theoretical calculations suggest that magnetic 
recordings can provide new information about the 
electrical anisotropy of cardiac tissue [3,4,5] and the 
problem only exists to reveal and use in clinical 
practice. 

The uniform anisotropy is due to two properties, 
which are microscopic in nature; the first being the 
difference in electrical conductivity in the axial and 
transverse direction of the fibers. Such anisotropy is 
described by the bidomain theory and characterized 
by the common fiber orientation defined by global 
angles in standard spherical co-ordinates [6], This 
theory explains experimental potential distributions 
measured from isolated canine hearts [7] and it was 
found that the conduction velocity parallel to the 
fiber axis was greater (3:1) than that of the 
transverse. Secondly, uniform anisotropy is caused 
by directional differences in cell-to-cell coupling on 
a macroscopic level. This anisotropy is described 
within the framework of the models of 
discontinuous conduction [8], 


Today ECG plays an important role in clinical 
diagnostics and cardiac research, and an equivalent 
electric dipole is used as the model of the cardiac 
electric source. The said dipole, in fact, is the vector 
of the electro-motive force (EMF), i.e. electric 
potential, which is also called the electric heart 
vector (EHV). 

The most used dipole within the MCG-method 
model is an equivalent magnetic dipole which is also 
called as magnetic heart vector (MHV) and 
equivalent current dipole (ECD) which is the vector 
of the impressed current. It was shown that angle 
between EMF and MHV at the QRS peak is very 
close to 90° [9], It should be mentioned that if this 
angle is exactly equal to 90°, there would be no new 
information in the MCG compared to ECG. 
However, it has been found that this angle varies 
considerably during QRS complex, both from 
patient to patient [10] and within various cardiac 
disorders [11]. 

The above-mentioned results show that the angle 
between MHV and EHV should contain new 
information about the myocardium conductivity. 
However, in our opinion MHV does not reflect 
completely electrophysiological phenomena because 
the physiological sense only has the multipole 
expansion of the cardiac generator [12], It is sum of 
the equivalent current generators-multipoles where 
ECD is the equivalent current generator of the 
lowest order. 

In the overall anisotropic case a myocardium 
conductivity is a tensor [13], However, at the axial 
symmetry with respect to fiber direction one can 
consider only two components - transversal (Sj_) and 
axial (So) conductivity. The present study is based 
on the statement that the degree of anisotropic 
conductivity manifestation is characterized by an 
angle y between a direction (a) of a EMF vector and 
direction (P) of maximal current density (MCD) 
vector in the human heart at the same time 

Y=a-p. (1) 

Then, isotropic I 0 and anisotropic I_t components of 
ECD vector J take the form 



fo = J cos y 1± = J sin y, 
tany=S±/S 0 . (2) 

Within the bidomain theory [3] one can identify the 
I o as isotropic current dipole and Ij. as that of 
anisotropy. 

The aim of this study was to investigate the 
connection between the degree of anisotropic 
conductivity manifestation (ACM) and myocardial 
ischaemia in patients with coronary artery disease at 
normal resting ECG. 

2 Methods and Results 

We investigated 27 patients, with CAD (stenosis 
>75% in > 1 vessel), stable angina, without 
myocardial infarction and with normal resting ECG 
or unspecific changes of ventricular repolarisation. 
The control group consisted of 18 healthy 
volunteers. 

The direction a of a vector EMF was calculated at 
the R- and T-peaks by standard ECG data. 

The MCG was registered at a rectangle grid of 6x6 
with 4 cm pitch by means of the single-channel 
magnetocardiographic system (MCG1, SQUID AG, 
Germany), installed in an unshielded location. The 
above device provides the registration of the vertical 
magnetic field gradient component d^B/dz^ and 
reference ECG. The said signals continuously 
sampled under the discretization frequency of 1000 
Hz over a period of 30 s at each point of a grid. 

The software provided sampling, filtration, artifact 
removing, averaging and magnetic field distribution 
maps calculations. 

The program CURLINES realizes the method of 
magnetic moments (MMM[14], According to this 
method a layer of magnetic dipoles distributed on 
the investigated surface is considered as heart 
biomagnetic field sources. In this case this surface is 
a plane located inside a thorax. Using the MMM- 
method, the correspondence between the unknown 
distribution of layer density in the sources plane and 
the distribution of a magnetic field induction of 
heart received in the measurement plane, is 
described by appropriate integral equations. The 
solution of these equations results in a distribution 
of a magnetic charges double layer density (or, in 
other words, function of a current), the level lines of 
which present lines of a current. 

The reconstruction of these currents is displayed in 
the forms of current density maps and current line 
maps. Current lines are curves along that the current 
flows, i.e. at any point the elementary current directs 
tangentially above line. The number (spatial density) 


of current lines is proportional to current density at 
that point. The spatial point of the maximum current 
density and direction of vector, which is tangential 
to current line at that point, are determined. 



a) 



Figure 1: Mean ±SE values of angles y: a) R peak, 
b) Tpeak. 

The vector reflects the direction of the maximum 
current density (MCD) used for calculation of the 13- 
angle at the R- and T-peaks. After calculation of the 
P- and a-angles, the value of ACM was obtained by 
expression (1). From Fig.l one can see that at the R 
peak the angle y lies within the range 20.9 ± 2.82 
(mean healthy volunteers, and within the range 
47.52 ± 5.45 for CAD patients. Within the T-wave, 
the angle lies within ranges 19.25 ± 3.84 and 
49.8 ±7.2 for the healthy group and CAD pts, 
respectively. Thus, the above results suggest that 
angle of ACM increases in CAD patients in 
comparison to healthy volunteers. The above 








increase is slightly more for T peaks(by almost 2.38 
times) compared to the R peak (almost 2.27 times). 
From the Figure it is clearly seen that the ACM 
angles for both groups are significantly different. T- 
test analyses show both groups are significantly 
different; p<0.025 for the R peak and p<0.01 for the 
T-wave. Sensitivity of this index is 93%, specificity 
85%. 

3 Discussion 

The effects in the presence of both spiral and 
uniform anisotropies will be summarized. In the 
simplest case it leads to simple summation of 
respective partial angles of ACM (1). The healthy 
person is characterized by the summary angle 
(between the EMF and MCD) which are determined 
by both the spiral structure of the heart and the 
conductivity anisotropy. Spiral anisotropy changes 
from one human to another and depends only on 
anatomical peculiarities. The conductivity changes 
in a healthy person but it also changes due to heart 
disorders. The ECG method is insensitive to spiral 
anisotropy and weakly sensitive to uniform 
anisotropy, but the MCG method is sensitive to both 
kinds of myocardium anisotropy. Results [13,15, 16] 
suggest that the myocardium anisotropy cannot only 
be recognized by ECG. 

The presence of pathological variations in part from 
all numbers of myocytes, leads to a change in the 
average longitudinal component So of the 
myocardium conductivity. Therefore, disturbances 
of the excitation wave propagation will be 
characterized by changing the ACM angle. Thus, 
increasing the ACM angle (Fig.l) is explained by 
decreasing of S<? due to existence of lower 
conductivity zones under CAD. 

As propagation of the above wavefront has two 
qualitatively different phases, de- and repolarization, 
the comparison between normal and pts groups at R- 
and T-peaks was logically performed. However, the 
use of R and T-peaks are not restricted in the present 
approach. The calculation of the ACM angle can be 
performed for any cardiocycle timepoint. The 
difference between the ACM angle at R- and T- 
peaks could be explained based on the difference of 
the excitation front propagation in the course of 
depolarization (mainly radial) and repolarisation 
(mainly tangential). The second reason may be that 
the above phases are characterized by different types 
of ions mainly carrying the transmembrane currents 
(i.e. Ca ions for de- and K ones for repolarization). 
In this case, ACM may be connected with different 
influences of certain heart disease to conductivity 
ions of selected species (e.g., due to their mobility or 


membrane permeability). Thus, decreasing of the 
ACM angle at the T-peak in comparison to the R- 
peak may by explained by increasing the 
longitudinal conductivity averaged over 
myocardium involving repolarization at that 
timepoint. On the other hand, the same effect under 
CAD is also explained by increasing of the K 
conductivity or decreasing that of Ca. 

4 Conclusion 

The proposed approach significantly distinguishes 
the group of healthy subjects and group of patients 
with CAD. Thus, analysis based on determination of 
the ACM using both ECG and MCG data in 
combination, is potentially useful for identifying 
CAD with high sensitivity and specificity. 

It should be mentioned that the presented approach 
is oriented towards clinics equipped with very 
simple ECG and MCG apparatus. 
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